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The Diurnal Variation of the Atmospheric Potential Gradient 
on the Summit of Mt. Fuji and along Its Slope. 


By H. HATAKEYAMA & K. UCHIKAWA 


Meteorological Research Institute, Tokyo, and Honjo Branch 


of the Aerological Observatory. 


Abstract 


The atmospheric potential gradient was observed on the summit 
of Mt. Fuji (3778m), at 5.5 gd (2800m) and Gotemba (460m), each during 
two or three days in August, 1947, and January, 1948. For this 
purpose, a polonium-collector was attached to one side of Shimizu’s 
unifilar electrometer opposite to the side the microseope is fixed. 
This instrument was made by the Institute of Physical and Chemicai 
Research of Japan, and it was so designed that it can be mounted 
on a tripodv The visual observation was made hourly or half-hourly 
to investigate the diurnal variation of the atmospheric potential gradi- 
ent. According to the observational results, the ordinary semi-diurnal 
variation (i.e. two maxima and two minima) prevails at Gotemba but 
at Tar6dbo the diurnal variation is not so conspicuous, whereas at 5.5 
go and at the summit the potential is greater in the daytime and 
smaller during the night. Thus it was made clear that the large semi- 
diurnal variation of the atmospheric potential gradient shows a change * 
in character also in the vertical direction besides the local influences 
due to surrounding conditions. 
1. Introduction. 
It is well known among meteorologists that the atmospheric potential gradient 
observed near the earth’s surface shows its maximum value in the morning and in the 
evening. As to the cause of the phenomenon, Whipple’s theory” that the large value 


of potential gradient is in connection with the turbidity of the air near the earth’s 


surface is the most generally accepted. While for the upper atmosphere there is no 


- observation of the diurnal variation of the potential gradient available, except some 
- made on aeroplane®® and also by Lautner“ who made continuous observation, though 
_ for a short period, of the potential gradient on the summit of Zugspitze (2312 m). 


Investigations, in order to separate the world-wide variation from the diurnal 


_ yariation of the potential gradient near the earth’s surface, have been and are still 
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being done by atmospheric electricians. So it is very interesting to. know how the 
type of the diurnal variation of the potential gradient changes with height. The 
authors of the present paper observed the diurnal change of the atmospheric potential 
gradient with a portable electrometer on the summit of Mt. Fuji (3778 m), at 50°20 
(2800 m), Tardb6 (1290m) and Gotemba (460m), each during two or three days in 
August, 1947, and January, 1948, with Mr. S. Ohta, who made observations of the 


condensation nuclei. 


2. Method of Observation. 


The electrometer used in this observation was Shimizu’s unifilar electrometer 
made by the Institute of Physical and Chemical Research. The metallic case of the 
electrometer and the microscope were removed from the stand and mounted on the top 
of the tripod of photograyhic camera 
as photographically shown in Fig. 1., in 
which A is the poionium-collector and 
B is the ambroid. The metallic case of 
the electrometer, the microscope, dry 
cells, wires are put into a small portable 
wooden hox shown in the lower part of 


Fig. 2 with other accessories. 


As shown in ig. 1, the case of 
the electrometer is earthed with a thin 
wire and its height is about 1.2m from 
the earth’s surface. The polonium-col- 
lector is about 3cm apart from the 
earthed case of the electrometer. When 
the observation is made on the plane 
earth’s surface, the atmospheric poten- 


Fig. 1 Portable Electrometer A: Collector, tial gradient is obtained from the indi- 
B: Ambroid cation of the electrometer-by multiply- 
ing the ordinary plane reduction coefficent. An example of the coefficient is shown 
in Table 1, which was measured at Kashiwa.race-course, Chiba Prefecture, on September, 
1947. pat 
The scenes of actual e Table 1 


~ 


observation are shown in Height of the| Coefficient of the 


Figs. 2 and 3. Fig. 2 shows collector (m) | plane reduction Remarks __ 

that at Tarob6 and Fig. 3 EF — +2102 | For Tarobd and Gotemba’ 

that on the summit of Mt. 0.58 4.91 | For 5.5 go wits: “S 

Fuji. On the summit the 234 7.26 | Collector on the upper. corner — 
; 5 oem SORE | 5 ‘OL CA corn 

porentiiweridinintadetay 0 9.77 Collector oA. the lower. « er 


_large that we can hardly catch the image of the fiber in the field of microscope. 
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Fig. 2 At Tarobo Fig. 3 On the summit of Mt. Fuji. 


In the actual observation, there are growing trees in the neighbourhood of the 
observing spot as seen in Fig. 2, or the spot itself is a peaked summit as seen in Fig. 
3. Thus another reductiou factor to reduce these topographical effects is necessary. 
To determine the reduction factor we have no other means than to estimate it from 
the position of the electrometer and the topography around it considering the equipo- 
tential surface. But we do not apply this reduction to our observations because our 
aim is only to clarify the type of the diurnal variation of the potential gradient. 

For the observation during the night, electric torch was placed on the same side 
as the collector. The indication of the electrometer becomes smaiier when the electric 
torch is brought near by. The reading of the electrometer during the night was reduced 
to that of the daytime by multiplying the factor duly determined. 

In August, 1947, the observation was made at Tar6obo from 14h, 18th to 11h, 
19th, at 5.5. gd from 7h, 20th to 8h, 21st, on the summit from 16h, 2lst to 7h, 24th, 
at 5.5 g6 again from 12h, 24th to 9h, 25th and at Tardbd again from 14h 20m, 25th to 
9h 30m, 26th. The observation was made hourly in principle, and -occasionally half- 
hourly or bi-hourly. The oper reon was discontinued for several hours at night for 
sleeping. 

In December, 1947, and January, 1948, the observation was made at ES from 
lib, 29th to 11h, 30th, on the summit from 17h, 2nd to 11h, 4th, at Tardbd from 9h, 
6th to 12h, 7th and at Gotemba again from 15h, 7th to 17h, Sth. <.<=2 
3. Results of Observation. : 
| _ Results of observation in Summer is shown graphically: in ‘Fig. -4 and results in 
Winter in Fig. 5. The abscissa is the time and the ordinate- is the potential. gradient 

expressed i in volts per meter. As already mentioned the value of the potential gradient 
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to the pole-effect. 
As to the type of the diurnal variation of the potential gradient, that at Gotemba 
exhibits ordinary semi-diurnal change (i.e. two maxima and two minima). One maxi- 
munt is at 8h or 9h and the other is at 17h or 18h. From parallel observations of 


_ other meteorological elements, it was found that these maxima have close relation 


with the trailing haze-layer or the coming and going of the turbid air caused by the 


sudden change of the wind direction. 


At Tar6dbé the potential gradient is almost constant through the whole day both 


- in Summer and in Winter. 


At 5.5 go and the summit the potential gradient is greater in the daytime than 


. in the night, being almost twice as large. The type of diurnal variation at Zugspitze 


already referred to is similar to those at 5.5 gd and the summit. This type of diurnal 


variation is considered as universal on a mountain higher than 2000m. 


WAY) ert petnr 


Mr. S. Ohta observed the number of the condensation nuclei at the same time 


ke B with our observation of the potential gradient. According to his conclusion the number 


of nuclei becomes smaller as the height increases and they are most numerous just 


after the noon and least in the early morning. Also he concluded that this diurnal 


variation of the number of nuclei may be due to the transportation of the nuclei by 


the ascending and descending currents along the slope of the mountain. 


The potential gradient becomes greater (smaller) when there comes an air current 


_ containing numerous (less) condensation nuclei, consequently less (more) smail ions 


and having less (hetter) electrical conductivity. These are agreeable observational — 
results. In the free atmosphere | the range of the diurnal variation of. the potential 


gradient ‘will be less than that at the summit or ‘at the si de of a mountain, because 
. in the free atmosphere, the range of the diurnal variation of the number of nuclei . 
is considered less than that on a mountain, as Mr. + Ohta estonated Thee pone sal | 


be verified in the future. — = 
Sr 1 conclusion the aitanys. wish to “express their. warmest thanks 
Chief of M Fuji Weather Station, Mr. Ke Takei, Chief 


Sudden Changes in the Atmospheric Electric Phenomena 
accompanying Lightning Discharges (I )* 


By Minoru KAWANO 


Geophysical Institute, Tokyo University 


Abstract 


Changes in the atmospheric electric conductivity and in the space 
charge accompanying lightning discharge were examined during summer 
‘in 1948-1950. It was found out that the conductivity increases re- 
markably at the time of lightning discharge. We might be able to 
say that this increase in the conductivity is attributable to an anomalous 
ionization by some mechanism accompanying Ene. discharges. 
iE Introduction . Bee = 


It has been found out by several investigators (e.g. Learid-2 in References) that 


the marked increase in the atmospheric electric conductivity follows thunder-showers. 


Changes in the electric field accompanying lightning discharges were discussed in 
- detail by Y. Tamura (3). Those in the atmospheric electric conductivity and. in the 
space charge in the same cases, however, have scarcely been observed, In our ijabora- = 
== tory ees observations of —— in the spores soos the Sager hates = 


“changes in” the: atmospheric electric conductivity, aaa in Ane space nested wil be 4 ia Te 3 


nas : so ‘ 
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Fig. 2 


connected to the quadrant electrometer. The height of the collector is 1.0 metre from 
the earth’s surface. The linear velocity of rotation of the recording drum is 10 mm/min. 
Ill. -Results of Observation 

The lightning dischages were detected by the records of the atmospheric electric 
field measured by the mechanical collector. The delay time of this collector is about 
0.2 second, and the record obtained by it can be thought as well following to sudden 
changes in the atmospheric electric field. By this reason the time of occurence of 
lightning ieee es are determined considerably acculately. 


Digew 


ae 


oe tenn 


Aug, 24, 1950 
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all observations. The average .of the anomalous increses in the conductivity was — 


about 40-50% of that at the time of thunder-showers. 

On the other hand, the space charge changes remarkably with the approach of 
thunder-cloulds. Each change, however, does not always correspond to lightning 
discharge. 

IV. Discussion 

It was ascertained by the observations of the electric field that the thunder- 
clouds which give rise to the above-mentioned anomalous changes were always those 
of ‘positive polarization.’ 

Judging from the observed facts that change in the space charge accompaning 
lighting discharges were very small, the positive partial conductivity would probably 
increase simultaneously with the increase of the negative conductivity. If sudden 
changes above-mentioned were caused in the electric field at the time of lighting 


discharge, the space charge will also have to be changed, because the electric charge 


induced in the earth’s surface should also be subject to change. The sudden change 

in the atmospheric electric conductivity accompanying lighting discharge, therefore, 

seems not to be due to tee change in the electric field. 
To explain this phenomenon we must consider some Beataione ionization near 


— 


ee ee ae eee ee ae 


i A al al be 


the ground at the time of lightning discharge. If a high energy radiation from < a 


- lightning discharge is assumed, we shail be able to Sola this phenomenon reasonably. 


= E. (@ Halliday 4), in fact, did observe such high energy ~ radiation: SCCOnI PSE a aS 


ight discharge. 


ae may be -poted that the detailed examination of thie: phenomenon will he = 


a Ad cnowledgment re Rete ; : ras 2 


a Se ase the writer wishes to pele his sincere e thanks to Dr. is oe 
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Current System for Sp-Field and the Bay Disturbance* 


By Naoshi FUKUSHIMA 


Geophys. Inst., Tokyo University 


Abstract 


The dynamo-theory for Sq has been proved to be applicable also 
for Sy» or the bay disturbance in the geomagnetic field. The distur- 
bance field of Sp or bay is generally larger in the dark hemisphere 
compared with that in the sunlit hemisphere. This fact can be inter- 
preted as a result of the condition that the electric conductivity of 
the auroral zones is markedly large, especially in the dark hemisphere. 
In this paper, the current system derived from the dynamo-theory 
with such conditions is numericaily computed for various probable 
cases with the aid of the perturbation method. The calculated current 
system, such as given in Fig. 5, is in fairly good agreement with the 


observed one. 
§ 1. Introduction 


It has been shown by T. Rikitake? and further extended by the present author 
that the current system corresponding to Sp-field and the bay disturbance can be 
derived from the dynamo-theory under the assumption that the conductivity above the 
so-called auroral zones-is much higher compared with that in the other regions. In the 
calculation hitherto made, however, the conductivity in each region was assumed to 
be invariant with longitude as well as latitude for the sake of mathematical simplicity. 
Although the calculated current system well agrees, as a first approximation, with 
that derived from the actual analysis of the geomagnetic data, there remain some 
discrepancies between them; For example, the current intensity is large in the dark 
hemisphere compared with that in the sunlit hemisphere. Such a discrepancy might 
be removed by adopting the more plausible assumption on the distribution of the 
conductivity of the upper atmosphere. In this paper, the calculated current system 
under the appropriate assumption of the conductivity was compared with the observed 
result. The method used here is a similar one with that used by A. Schuster® in the 
case of examining Sq current system, that is, the perturbation method. ==" 


* Contribution from the Division of Geomagnetism and Geoelectricity, . 
_ Geophysical Institute, Tokyo University. Series II, No. 12 (1950) 
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§ 2. Assumptions on the distribution of 
the conductivity in the upper atmo- 
sphere. 


As shown in Fig. 1, the conducting 
layer of radius Ry, which is concentric 
with the earth, is divided into 5 parts 
with respect to geomagnetic latitude, 
namely 


the polar caps (regions I and V) 


m= SS IV (@=0°~20°, 160~180°), 
i LLLLILLL the auroral zones (regions II and IV) 
(@=20°~25°, 155°~160°), 
Fig. 1 and the equatorial zone (region II) 
(@=25°~ 155°). ‘ 
In the equatorial zone, the total po sca ane Sn of the conducting layer. K is 


assumed to be 
K=K, (1+7’cos x), 
; where A, is a constant, x denotes the zenith distance of the sun, and nee constant _ 
Tess than unity. If the sun’s declination is zero, 
= cos x= —sin @ cos ¢ 
where ¢ is longitude reckoned from the midnight meridian. 
As for the conductivity in the higher latitude regions, its detieadentes on x oonita 
a “be ignored, because 7’ cos x€1. According to the result of the actual analysis of ‘eal 
ee. "geomagnetic data in the auroral zones, such as by M. Hasegawa®, L. Harang”, HC. | 
__ Silsbee and EH. Vestine, the conductivity there can be considered to be very large 


ae penny. in = ae hemisphere. It will Lbs recognized that its maximum value takes 
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¢, and ¢ is assumed to be expanded by surface spherical harmonics, namely 


g= WHkyPx (cos 6) sin (My+a%). (2) 


The earth’s magnetic field is replaced by the magnetic field of a centred dipole of 
moment M, with its axis coincident with the rotation axis of the earth. Vertical 
component of the earth’s magnetic field Hz can be given by 


region given by the relation (1). 
3. 1. Solution in the auroral zones 

Substituting the expression for the conductivity in the auroral zones into (4), and 
negiecting the higher order terms of 7, (4) becomes 


Hz=— 8 cos 6=—2G cos 6. (3) 
Ro 

Then, the current function J for the stationary state must satisfy the following 

_ equations: 
fae 0s 

x Ro sin 6 0g aie ER a0” (4) 
Ss ay as 
<i -- = —KuH;—k —————,, 
t R06 ues R sin 60¢ 
b where S denotes the electrostatic potential, and K is the total conductivity in each 


| 4 | 2bKoG cot @ —— =bKy ae (1—y cos ¢) — ey 5 
sinf0¢g ; 
| ope (5) 
~2K,G cos 20 bK, —*S__(1—reos)- 2E., 


“0 06 sinfd¢g 
Tet So and Th denote the expressions for S and J when BS =o respectively. Then, in 
general expression, 2 
| geteees ie Esty ae : s'=0() } 
eer peop eto 
Where S and J satisfy the following equations : 


eS he Ay 


Foes ee 46) 


oe Ie SELOSY ry hes ; : ape 
eee ies Ee See ae 
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where 
4 uy N {(N- t1p— Me} (N+1)(N?—M)? rh 
S(O) =k’ a (V+1(2N+1) P*,,(cos 0)+ NON) P¥_(cos 0)| (11) 
Pn (0)= sin" 14+ Nes )) cos” 6 a. vee 


rrr oe, -{(2S—2)(2S— 142m) + m+ Dhogeesg 4... | ae 


(2S)! 
Qm(9)=sin”0| cos@ + otemt ant) costo +: -: 
4 2+2m-+m(m+1)}{34+2m)+m(m+D)}--{CS—VES+2m)+m(M+V} ossvag 4 .. Ja 
@S+)! | 
and C.”, Cy”, Co”, Cr", aft, Gi, Gu and af zy are the arbitrary constants, The | 


first term in the right hand side of (10) is the particular solution, while the remaining 

two terms are the general solutions. The particular solution of J; should be substitut- 

ed into (8), and then J’ can be obtained. 

3. 2. Solution in the equatorial region 
Substitution of the expression for the conductivity in this region into (4) gives — 


the equations : 


as eta 2 - Se 


2K Goor Se =K, 95 +(1 47’sind cose) ~- : 
Be = = 0S = Fes ao ; = | ; 
2KoGcosé 30 =Ky einer (1+7’sin#cos¢) E oS = 2 


which pemepand to (5) in the above case. After a ee treatment as before, the 
ee for J in this region will be written as follows : = Eee. 


= biG Sey Se ne = 
90 0g sinddg — 3 big Bee a 
as) ~faintse Mtoe sa ab 


pane: 


107 


and the solution is 
j= —2aKoG Sy (0) sin (Mg+a%)— 2aKoG dCi" (0) sin (my+a;?”) 
for the northern cap, 


fe. ~2aK GD) Sy (0) sin(M¢+a¥)— 2AK0GYCH"Un (0) sin (me+a™ Ge 
for the southern cap, 
where 
— ein ea m (m+ 1) ES oe 
Um (0) = sin” 0 E 2-1! G41) tnt) (cos 0—1)+ 
sm (m+1){2+2m+m(m-+1)} +: {(S—1) (S+2m)+m(m-+1)} “4\ 6 
+(-)) ane (cos 0 —1)5-+ | (20) 
Vm (0) =Um (x— 6) : ; ; (21) 


and C,”, Cs", af" and ay” are the constants. 


4 


_ 3. 4. Boundary Condition 

i The solution of the differential equation for the current function J is completely 

y © obtained in the above treatment. The groups of arbitrary constants C;”, C.,...., Cs” 
and phase angles in the expression for J must be determined by the conditions that 

. _ the normal component of the current and the tangential component of the electric 

“field serve respectively by 


ei - gc Of-= 
Ry sin 609, K R,00 


are continuous at the four boundaries of the five regions. 
3.5. Separation of Sq field and Sp field 
I will be convenient to consider that the Esper for the carrent function Se 
consists of two parts, that is, ge 
$ See J=Jq+In a i et ee Se 2) 
; "where the current system deduced from Jq corresponds to 0 Sy field, while that from : 
Jo causes See Parther,, if, # Jao and Joo denote the Ses for Jga and = when 2 = = 


ernie vectey | 
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in region I J=—2aKoG ki eo PAD. 4) = + Cy Uy (o)| sin (g+a1') 


+ {7Di+7'Dy } ts (0) sin (2 g+o)+{7 Bit7' By }sin a | 


a b] 
eeefecion Il. aif = = abkele bee 


a {(a (+ 5 & 9 pe () + sin 20)+ (7De+1/D'») pO) 
+(7D3+7'D's) qe (0) |sin2y-+os" 


4 fy Boia 04 (Be-tr!B')bo (+ (B+1'Bs) a0) | sin at | 


in region II J=—2KoG hs 


+ Cops (0) + Coa @| sin (¢+o1") 


{ Pe! (cos 6) 
L 2v3 
ae es ie ; -sin? @ cos 045% O))+(D. 28 D's) po(0)-+(7Ds-tr! D's qo) 
; X sin (29 +411) 
+ fare 30—3 cos 0)+(7Bs+7'B's) po ie + (Bs +7'B's)qo (6) \sin a’ | 


+ Cubs (OV Coat @| sin (g-+a1") 


\e5) 


i rebion TV fe 0K rs | £ a P3608 0) 4 cy, (0) + Cras (@)\ sin (y +01") 


— 


+ f(a +3 (F-0 )po(0)-+Esin20)-GDitr Do be) 


= . +(Dr+7'D! 7) ()| sin (2 +a) 


be = (gain 4 oer rFOMO-Bery Boa an a) 


Geran Cs v1 Ojsinte (y +a1') 


ae = = cas 


= + frp Det Ds ve (sin sae 
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in the polar cap, 


1 
—2a kG ky ae ep a —+ Cy 01 (0)| sin (y+ax') 
in the 


b—1 
me) =28 KG bat * , £ ED serosa (0)| sin (o-+ a1" auroral 
zone, 


—2Ko G RC; qu (0) in the equatorial region, 


—2akoG ki! iD ats (0) sin (29+ a,')+ By sin ai'} 


~20K,G kt | {a OG —0 )pu (0) +sin 2.0-+ Ds po (0) +Ds qo (0) 


X sin (2 2 obo!) +1 Lin 20+ B.—B; log tan.) lain a 4 in the auro- 
(8 25 ral zone, 


in the polar cap, 


in the equatorial (28) 


—2 KoG ki! | Dua (A) sin (2¢9+a1')—B; log oo sin a1 f . 
2 region, 


’ —2a K)G al {Dp ' Ue (0)—57- (cos 30—cos O)\sin (2 p+as') + Bi = 


x (cos 30 ed sin as*| in the polar cap, 


Toy \ —2b KG asf [ [D+ Ds) aay (008 30—cos )\sin (2o-+a1" 
in the auro- 


So ‘2 2— B’; log See 3450S 30— ace so) sin at | 
2 J ral zone, 


| ies cece elie ee 


— D-*= 2(0) sin 2¢+ a1')— B's log tans-sin as*| in the 
= equatorial Tegion, 


= 


in the eres hemispheze. Cy D. Ds a Be and By in (25) ar ce all zero. 
—=—— ‘In the he numerical calculation here, 


| 5 or 
t= : 2 es Se K, is assumed to be 3x 10-* em.u., M . 
~ 79x10" cgs., Re 647 x10" ¢ cm, that ise 


the ronda 
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intense in the dark hemisphere and less in the sunlit hemisphere, while the equatorial 
current intensity shows only a slight decrease in the dark hemisphere and increase 
in the sunlit hemisphere. The ratios of the total amount of the current flowing in 
the auroral zone and the equatorial region of the dark hemisphere to those of the 
sunlit hemisphere, denoted by f(A) and /(E) respectively, are given in Figs. 3 and 4. 
At that time, the current intensity in the polar caps does not show remarkable change, 
its change being less than 10%. 


01 02 03 04 05 06 y 01 02 03 04 05 06 » 
Fig. 3 Change in the ratio of the ~ ‘Fig. 4 Change in the ratio of the _ 
intensity of the auroral zone ‘current intensity of the equa- 
current in the dark hemis- : torial region in the dark 
phere to that in the sunlit aoe hemisphere to that in the 


hemisphere with y. sunlit side with y. 
In Fig. 5, the current system when ree and 7/=0 is shown. The current system for gz 


the bay disturbance obtained by H.C. Silsbee and E.H. Vestine” is also shown here for _ 


— 
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30° and of the current intensity between them. The latter difficulty may be removed 
in some way, for example, by the author’s hypothesis of the double conducting layers. 

In the analysis of the geomagnetic disturbance field at the time of magnetic 
storms, the Dst-component should contain, according to its definition, not only the 
magnetic field which is attributed to the equatorial ring current, but also the zonal 
part such as derived from J’py. In other words, the zonal part is to be eliminated 
from the so-called Sp-field. If this zonal part is eliminated from the current system 
shown in Fig. 5, the current system shown in Fig. 7 can be obtained. In. the latter 


system, the current intensity in the sunlit hemisphere and that in the dark hemisphere 
are nearly the same. 


ee a cae 
Fig. 7 The current system obtained by 
eliminating the zonal part from the 
current system in Fig. 5. _ 


hemisphere and 


112 


more complicated, but look more like the observed system. ; 
It must be noted here that the conductivity of the equatorial region Ko, which 
is assumed conventionally here 3x 10-® e.m.u., should be replaced by the more reliable 
value of about 5X10-%-e.m.u. obtained by T. Nagata. However, we concern in this 
calculation under dynamo-theory only the product of the conductivity and the velocity — 
of ionosphere wind, absolute values of conductivity being not significant. In other 
words, the relative magnitudes of conductivity in respective regions are to be the 
cause of Sp- or bay-fields under the consideration of the dynamo-theory. 
; At any rate, the dynamo-action is considered to be the most fundamental process 
to produce the various kinds of geomagnetic variation fields, although there still 


remains some difficulty in detail of each variation fields. 


In concluding, the writer wishes to express his hearty thanks to Dr. T. Nagata 
for his direction throughout this study, and to Miss N. Ono for assistance in the 


numerical computation... 
(Read: Oct. 9, 1950) 
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On the Influence of the Hall Current to the Electrical 
Conductivity of the lonosphere, II 


By Motokazu HIRONO 


Geophysical Institute, Kyoto University 


Abstract 


It is shown that there is a narrow region along line of zero dip 
near the E layer, of about 15° of latitude wide in which conductivity 
is very great, as a result of polarization by the Hail current. The 

effect of this belt is discussed by the dynamo theory, giving the result 
that this belt will cause enhanced diurnal variation of magnetic field 
near the geomagnetic equator. 


: Introduction. In the previous paper [1], hereafter referred to as I, we made 
a preliminary calculation on the influence of the Hail Current (in place of wich we 


write H.C. in the following) to the electrical conductivity of the longsphert Fhe 


present paper describes a further examination of the problem. : 
TAG. Cowling [2] showed that, when iste are few negative ions in the sence: : 
electrons s eparately. Therfore, mean velocity of | 
eeparecd particles i is expressed in the following eaiation = 


th "ps C= ve. ex ws 


= as : or 5 


E=BO+GxH 5 hn 
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Cowling [2] called o’ and o”, direct and transverse conductivity respectively. In the 

following, as in I, we treat the case 4<1 near the E region. The ratio o%/s' with — | 

gas density 2 (particles/c.c.) is shown in Table 1. 

Table 1 When o”€o!, as in the upper part of the Fl layer and in 4 

: H/o! the F2 layer (n<10"), current flows almost in the direc- 

‘tion of E. But when o”>o!, as near the E layer (7=6X ; 

a fae 10"), then an impressed electromotive force FE, produces not : 

102 3.0 only direct current, but also a considerable H.C., and we 

1013 17.6 must consider the effect of the polarization field H; by the 3 

10" 4.1 latter current. Then H=H,+41, provided # is perpendi- 

cular to H, therefore equation (5) becomes 
J=0'(Ey)+E,)+0Uh x (Eo +E) (7) 


rei ie > 


2. Conductivity of the EF region. : 
Let right-handed rectangular axis o-xyz be taken, so placed that ox, oy and oz 
-are directed to the south, east and upward respectively and unit vectors in the direction 
of ox, oy and oz be #, 7 and k respectively. According to the “ Dynamo theory ” 
concerning the diurnal variation of the geomagnetic field, both horizontal air velocity 
¢) and magnetic field A, produce primarily dynamo field ce x H=vH. t—uH.j—vHk. 
=e Schuster [3] in 1908 showed the possibility that —vH,k night produce a three 


_ dimensional current system. But we think that, as the thickness of the ionosphere, 


= “even when E and F layers are combined, is sufficiently thin compared with the kori- 
~ zontal uniformity of the diurnal current system, —vH,k is counteracted by the resulting x E 
electric field Ek and —vH,+E,~0. When we discuss the diurnal variation, therefore, 
we can eliminate both —vH,k and E;k from E in equation (1) and agree with Cowl. 
‘ing’s ‘opinion [4]. is I we discussed the effect. of the vertical component of the HC, 
in the oes, we discuss effect of the horizontal af component of that i in ssontoy: Asi in ts 


ae ctivity n near ihe feomagne- 


‘layer near 
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Substituting equation (8) in the upper equation, we have 
E,;=—(6"/o")-E) from this relation and. (8) it follows that 


J=cl. EF, where o“=¢!{1+(6"%/a?)?}. (10) 


o is considered to be resultant direct conductivity, and identical with equ. (11) in I 
for ¢=0. 
In Figs. 2 and 3, o!/n. and o/o! are shown with 2 and for 4<1. 
" 
10 to! 
, t 
12 
q . 10 107 
es 10° 10” 
| aay eee ! q : 
; SAG 8ex\0" = 10. Soe eee 
— oF /Ne (emu. bit : 
Sere a2 Fig. 3 3 
From these figures it is evident that o/”/n, is maximum in the lower “part of the E- 


(ii) Conduc tivity in middle latitudes. — Se = — = 


— 
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polarization electric field resulting from H.C. be #; and as current between AD, we 


have, using equation (1) : 
J=0(Fy+ Bi)" +o"hx (Eo+E1)+0 Es" where E,=Ey,i+(Eut Ek 

In I, vertical componet of H.C. and the resulting E’;. were discussed and the effect of these 
was shown to be negligible for conductivity in higher latitudes than 10°, therefoce we ~ 
ignore them in the following. 

Using the above equation and k-J=0, we dininnta F,. and get 


J=(o'E\,/sin? 6+o%E, sin $) t+ (o'E)—o"E)) J (11) 
where £,~£;,/sin ¢ 
e. b- J=o'E;,/sin2?d +o" sin 6+ Eo (12) 


where ¢ is the dip of geomagnetic field. 
As current density I integrated with height in every layer, we have for AD approxi- 


_ mately I= \ ol dh: Ej,/sin?bm+ \ ,o%-dh-E, sin $n, where $m is the mean value of in 
AD. Considering that, if AD, BC are not much less than AB, CD, the resistances bet- 
ween AB and CD are far less than those between AD and BC, and that, in the F 
region, as o”€a', even when similar Ey is. impressed there, its influence to i-J is 
negligible, for BC we have 

| \ yo dh: Ey,/sin?dn 
From above two equations a RTT dm? Eo eit =e = Fah Lae 
where a=—\,o%dh/(\, o'dh+), ot-dh) i agi 

Substituting (13) in (11) and denoting mean value in the F layer ore suffix: B.S 3 

: Ss a+o%)sin bm Ett (94, —ola sin? bm)EoJ eee nt eet : 

fe Jv=o%E, where O7} = Op A°SiM? bm Oy chet a 

: As Saoneit of J in i flows in an opposite direction for E and F region, their effe 
to the magnetic field is negligible. 

ee | a<1, Mey 15X10%, aa oldh= AI =I 
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lines of force, which crosses E layer, 


will not enter into the F region as 


3 


shown in Fig. 1 therefore o, will not 
greatly differ from o”% But at lati- 
tudes L>7.°2, as BC of the above men- 
tioned element ADBC enter into the F 
region, oy, will decrease rapidly from 
L=7.°2 towards higher latitudes tend- 
ing to 207. 


= 


~ CO Up to the present we considered 
—>{, the effect of H.C. for an element AD- 
CB. According to the principle of 
Big.6 superposition, the above mentioned 
Result will hold for the diurnal current also. When an electromotive force Ey)=(Ko. 
Evy) is imposed, we can take I,=06,-E),, I,=oy+Eyy with respect to the electric current 
producing the diurnal magnetic variation by reason mentioned in (ii). At the equator, 
for the mean value in the EF layer, in which local scale height is 10km, Syn=a1 it and | 
‘using equation (9) o//=400'. As L increases, oy will slowly decrease and beyond L= 
7.°2 oyx will rapidly decrease tending to 20%. Under the assumption that vertical 
ponent of H.C. does not flow, and using equ. (11) in I, we get for (ay/o7), curve 


A in Fig. 7. Such an assumption will be fit, when the effective layer is sufficiently thin —— 


CO! npared with current system. Curve A is considered to be a fairly good approximation 
or the above consideration, but curve B Sotresponding to eee will be better 
} for the latitudes L>7.°2. ae | 
“¢ mean value in 0°<L<7. °5 of oy Sec ee $e Sa Sa 
om a and that in middle latitudes 40/\ +t fo), 


= 


be P Fa then for pee actile uniferm = 
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current function be R, electrostatic potential be S, then we have the next relation [5], 
SOS eer OR 

aoe K asin 00¢ 

OS aes nett Ole 
asin 00¢ K ado 
where a is the radius of the shell, K is conductivity integrated with height. We assume 
the next distribution of K, I K=K;, (const.) for 0<0<0, Il K=nK, for 0:<0<n—-h;: 
Il K=K, for 1-0:<059/2 
We take, as the velocity potential of the air, for mathematical convenience 

v=k-P;'- sing (16) 

‘and investigate the outline of the influence of region IL.- 


(15) 


—uH.— 


SE OE Oe ee ee ee ae ee re eS eer 


Eliminating S from (15) we have 


eR 8 aR OH) , 0 (uH. sin 0) ) 
sin 00g? OA” G0. = ak} ELLIE ue 
We write R=dydNR", m = Pe Sin Mo. 
A particular integral of equation (17) in regions I, IL is R2'=(1/6)CKik P2'sin ¢, in region 
Il is Re! =(1/6) Cukyk Pe! sin ¢ © = (18) 
SS where H,=C-cos0. Next we consider the solution of the aeaaon 2g 
Spe ee OtR Pes O ORY = = ; er 
= Sin bag? Soe a sind 00 =; | , Boe (19) 4 
. Let R=f (0)-F (¢), MP be constant of separation and Sa gives. = 
= OF /0¢?=— FM = | Ose 
Sere = Of ee ee 5 $5 oe 
ge : an os in j= a0 =. 3 : | ; ja ep aa 


sin Mé 


Tn order th t Fb 1 = 
| a e sing e valued function of ¢, M must ie integer aul Fis ~ COS SMe 
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From above equations we can determine eastward current V, 


OR . 
= FERRE V, (@)-sin ¢. V, (0) is represented in Fig, 8 in unit of (1/6) CK, k for 


various values of 6; and yp. 


ee, MeL 


Fig. 8 


The more y increases, the more total amount of current increases,,and the ‘more 
current focus shifts towards equator and the more the current in lower latitudes con- 
-centrates towards equator. But, as it was shown by D.F. ‘Martyn that ‘some - diurnal 

currents flow in Fl and F2 layers [7], it is impossible to determine , “ eee from this 
. es ees the value. 2s = aught ‘be approprte: : 
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LETTER TO EDITOR 


On the Analysis of Sudden Indrease of Cosmic Ray associated 
with Solar Flare 


The sudden increase of cosmic ray intensity associated with the solar flare is of 
very interest because it seems to give some clue to the origin of usual cosmic rays. 

The questions which arise from this phenomenaon may be divided in.the follow- 
ing three parts: 

(1) Did what kind of incoming particles give rise such unusual increase of cosmic 
rays ? 

(2) Is there any possible mechanism that can accelerate the charged particle up to 
several Bey at the time of solar flare ? 

(3) If such mechanism is possible, can such particles escape from the strong 
magnetic field of the sun ? 

Questions (2) and (3) have already been studied and answered affirmatively by Swann“ 
and Forbush et al® respectively. 

In order to study the feature and the energy distribution of the incoming rays 
which produced such increase, the analysis was made using the data of the increase 
on Feb. 28, and Mar. 7. 1942, and the following results were obtained. (The available 
data for this analysis were obtained at Huancayo (1°Sm), Christchurch (48°Sm), 
Cheltenham, (50°Nm) Godhavn (78°Nm)® and Tokyo (25°Nm)). 

(1) The increase at Godhayn was larger than that of Cheltenltenham. The 


_ difference was surely caused by the particles with energy less than the magnetic cut 


_ off energy at Cheltenham. Furthermore, the particles must have at least 2 Bev per 


Poe 


nucleon to penetrate through the atmosphere. From Table I the platicles satisfying 


Table. 1. Magnetic cut off energy per nucleon at Cheltenham. 


particle | electron | proton | a particle | other nuclei 


Magnetic cut off 


~ 4 BeV ~ 3 BeV ~ 1.2 BeV ~ 1BeV 


energy /nucleon 


above conditions are protons or electons. However, the possibility that incident parti- 
_ cles were only electrons can be excluded from the fact that the increase was also 
observed under the thick (10 cm or more) lead absorber“. We thus conclude that 


almost of alli ncident particles were protons. This conclusion was confirmed by recent 


; experiment of photographic plate carrid out by oe et al, who found that So! 
; Protans: were increased at the time of solar flare, 


- 
' 
' 
+ 


a 
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If the particles had been accelerated so gradually that the heavy particles had 
not been destroyed, . they would have been accelerated and observed. From _ this 
argument it may be concluded that the particles were not accelerated so gradually. 

(2) From the latitude effect of the increase, the energy spectra of incoming 
particles which produced such increase are roughly estimated.©. As shown in Fig. tS 

~ the average energy of particles de- 
ue creases with increasing time on Feb. 


410 - 44 28. ; 
Yay. Feb. 28, while it is almost constant during 


Le ee ee ee ee a 


the increase on Mar. 7. : 
This contradiction may be ex- 


plained by the following assumptions : 
{12-16 py, Feb 28 (1) The tunnels through the 
Stérmers forbidden region opened up 
during the acceleration on both days. . 
While the acceleration mechanism was 
(4-18, feb,28. gradually demped with time on Feb. 
28, but was suddenly damped on Mar. 


02 


al a = 
7 after having accelerated the particles 


for a short time. 


é Sar ee 5 Baas : 7% ze 
Boe, e ; eee _ (2) The acceleration mechanism 
area Fig. 1. Differential Spectra of Incoming _ ; 
“See Particles. <10—14>AV. means the - While the tunnel opened up during the © 7 


was gradually damped on both days. 


a i average value from. 10 to 14 G.M.T. 
ee g acceleration on Feb. 28, ee on Mar. — 


; assumption @. However, for oe lack of “information available to us, we. can aes 
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